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ABSTRACT 

Pristine stars with masses between ~140 and 260 Mq are theoretically predicted to 
die as pair-instability supernovae. These very massive progenitors could come from 
Pop III stars in the early universe. We model the light curves and spectra of pair- 
instability supernovae over a range of masses and envelope structures. At redshifts 
of reionization z ^ 6, we calculate the rates and detectability of pair-instability and 
core collapse supernovae, and show that with the James Webb Space Telescope^ it is 
possible to determine the contribution of Pop III and Pop II stars toward reionization 
by constraining the stellar initial mass function at that epoch using these supernovae. 
We also find the rates of Type la supernovae, and show that they are not rare during 
reionization, and can be used to probe the mass function at 4-8 M©. If the budget 
of ionizing photons was dominated by contributions from top-heavy Pop III stars, we 
predict that the bright end of the galaxy luminosity function will be contaminated by 
pair-instability supernovae. 

Key words: supernovae: general - dark ages, reionization, first stars - stars: Popu- 
lation III - stars: Population II 



1 INTRODUCTION 

The life of a massive star ends in a supernova (SN). The 
detection of neutrinos from SN 1987A verified the idea 
that some SNe are set off by the gravitational collapse of 
the iron core of their progenitor star (jKrausa.1987. ). How- 
ever, theory predicts that very massive stars with helium 
cores between ~64 and 133 Mq could find another way 
to blow up, through the thermonuc lear explosion of oxygen 
via the pair-production inst ability (Raka vv fc Shavivlll967l : 
iBarkat et a"lll967l : lHeger fc Wooslev 2002 ). The production 
of electron/positron pairs in the core softens the equation 
of state, leading to collapse and the ignition of explosive 
oxygen burning. The subsequent thermonuclear runaway re- 
verses the collapse and ejects the entire star, leaving no rem- 
nant behind. Unlike iron core-collapse supernovae (CCSNe), 
which involves poorly constrained physical processes such as 
turbulence, pulsations, perhaps rotation and magnetic fields, 
the physics involved in pair-instability supernovae (PISNe) 
is fairly wel l understood and can be modeled with fewer un- 
certainties (|Langeij|2009l ). 

Due to the extremely large stellar mass required, the 
progenitors of PISNe are expected to be rare, and may only 
form under unusual conditions. One such condition existed 
in the earl y universe, when metal-free Population HI stars 
were born ((Locb' 2010). In star formation, it is the accretion 
process that ultimately sets the final mass of a star. From di- 



mensional arguments, the mass growth rate is simply given 



by the Jeans mass Mj 



over the free-fall time 



^// ~ implying dM/dt oc c?s/G cx , where the 

sound speed Cs ~ \/ kT/mp. In present day star-forming 
regions, heavy elements radiatively cool the gas to a tem- 
perature as low as T ~ lOK. However, in primordial clouds, 
the primary coolant at low temperatures is molecular hy- 
drogen, which can only cool the gas to T ~ 200 — 300K, 
implying an accretion rate higher than present day by two 
orders of magnitude. Hence, theoretical studies suggest that 
the initial mass function (IMF) of Pop HI stars might 
have been biased toward masses much higher th an today, 
e.g. several hundred Mq l|Bromm fc Larson: '2004). The nu- 
cleosynthesis imprints of this top heavy IMF have been 
seen in globular cluste r s and damped L yman alpha sys- 
tems (|Cooke et al.ll201ll : |Puzia et al.ll2006l ). Moreover, mas- 
sive stars have strong winds driven by radiation pressure 
through spectral lines, with a mass loss rate scaling with 
stella r metallicity M oc ^O-S-o '^ (|Vink et al.l200ll : lKudrit^ 
l2002h . Most PISNe should therefore be from Pop HI stars, 
which have weak radiation-driven winds due to their ex- 
tremely low metallicities, and retain enough of their initial 
masses at the end of their lives to undergo a pair-instability 
explosion. 

Naturally, studies of the rates and detectability of 
PIS Ne focused on high redshifts before reionization. 
iMackev et all l|2003h found the PISNe rate to be ~ 50 
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deg"^ yr"^ at z > 15, while IWeinmann fc Lilly! (120051 ). us- 
ing more conservative assumptions for the number of PISNe 
produced per unit Pop III stellar mass formed, found the 
PISNe rate to be ~ 4 deg~^ yr~^ at similar redshifts. As- 
suming that only one supermassive Pop III star forms in 
une nriched minihalos, a nd that none form in protogalax- 
ies, IWise fc Abeil (|2005l ) found the PISNe rate be ~ 0.34 
deg~'^ yr~^ at z ^ 20. After our paper was submitted, 
iHummel et"al] (j201ll ) presented a complimentary analysis 
of the source density of PISNe from pristine minihalos, and 
determined the observability of such events with the James 
Webb Space Telescope (JWST), finding approximately ~ 0.4 
PISNe visible per JWST field of view at any given time. 
PISNe after the epoch o f reionization were also considered; 
IScannapieco et aD (|2005l ) calculated a suite of PISNe model 
light curves with blackbody spectra, and analyzed the de- 
tectability and rates of PISNe from Pop III stars formed 
from leftove r pristine gas at z < 6. As for CCSNe during 
reionization, iMesinger et al.l l|2006l ) presented detailed pre- 
dictions for the number of core collapse SNe that JWST 
could observe as a function of different survey parameters. 

In this paper, we present light curves and spectral 
time series for PISNe from our multi-wavelength radiation- 
hydrodynamics simulations. As the stellar population re- 
sponsible for reionization is currently unknown, instead of 
predicting a fixed SNe rate, we normalize the star forma- 
tion rate by requiring that enough ionizing photons must be 
produced by either Pop III or Pop II stars in protogalaxies 
to complete reionization by 2; ~ 6, and calculate the rates 
of pair-instability, core-collapse, and Type la SNe and their 
detectability with JWST for these two scenarios; the actual 
SNe rates will be in between these limiting cases. We show 
that using the observed rates of these SNe, it is possible 
to distinguish the contribution of Pop III and Pop II stars 
toward reionization by characterizing the IMF at that time. 



2 LIGHT CURVES AND SPECTRA 

The stellar evolution and explosion of PISN models, and the 
resulting broadband light curves and spectral time-series are 
described in detail in lKasen et al . (2011); here we summarize 
the resuhs. Models R150, R175, R200, R225, R250 represent 
explosions of 150-250 Mq red supergiant stars, respectively, 
each with their hydrogen envelope intact. In principle, blue 
supergiants are also possible progenitors of PISNe, but con- 
vective mixing of metals into the hydrogen envelope makes 
it more likely that the progenitor dies as a red supergiant. 
Models He80, HelOO, Hel30 were explosions of 80, 100, 130 
Mq bare helium cores. Such models may represent stars that 
lost they hydrogen envelope due to a prior pulsational phase 
or through binary interactions. Here we use an approximate 
empirical relation between the helium core mass and the 
progenitor main-sequence mass (|Heger fc Wooslevll2002l ): 

MHe~^(MzAMS -20Mq). (1) 

Properties of all presupernova stars and their explosions are 
given in Table [T] 

A few days after the explosion, hydrodynamical pro- 
cesses subside and the ejected material reaches a phase of 
nearly free expansion. The energy powering the subsequent 



Table 1. Parameters of supernova explosion models. The R-prefix 
models refer to red supergiant progenitor PISNe, and the He- 
prefix models refer to the exposed helium core PISNe. The proxy 
core-collapse SN model (CC), a model Type IIP supernova, is 
shown for comparison. _Ro is the presupernova radius. Mi, Mf 
are the initial and final masses of the progenitor, respectively, 
while Mhb is the helium core mass and Mffi is the amount of 
^^Ni synthesized in the explosion. All mass units are in Mq. 



Name 


M, 


Mf 


Mffe 


Mm 


Ro (lO^^cm) 


R150 


150 


142.9 


72.0 


0.07 


162 


R175 


175 


163.8 


84.4 


0.70 


174 


R200 


200 


181.1 


96.7 


5.09 


184 


R225 


225 


200.3 


103.5 


16.5 


333 


R250 


250 


236.3 


124.0 


37.86 


225 


He80 


80 


80 


80 


0.19 




HelOO 


100 


100 


100 


5.00 




Hel30 


130 


130 


130 


40.32 




CC 


15 


13.3 




0.28 


44 



light curve may derive from three possible sources: (i) Lin- 
gering thermal energy from the explosion itself; (ii) The 
radioactive decay of synthesized ^^Ni; (iii) The interaction 
of the ejecta with a dense circumstellar medium. Thermal 
energy suffers adiabatic losses on the expansion timescale 
tex ~ Ro/v, and so source (i) is only significant for stars 
with large initial radii Rq. Circumstellar interaction has not 
been included in the models discussed here. 

We have computed light curves and spectral time se- 
ries of the explosion models using the time depe ndent radia- 
tive transfer code SEDONA (|Kasen et al.ll2006l ). All models 
shown here assume spherical symmetry, and calculations of 
atomic level populations assume local thermodynamic equi- 
librium. Using Monte Carlo methods, we solve the full multi- 
wavelength radiative transfer problem using realistic opaci- 
ties as a function of wavelength, composition and tempera- 
ture, over mill ions of line transitions. Unlike previous black - 
body models (jWise fc Abell[2005l : IScannapieco et al ] l2005l ). 
our results allow us to calculate more accurate light curves 
for any given color bands and to study the time evolution of 
the supernova colors and spectral features. 

The shape and duration of PISN light curves depend on 
the mass and radius of their progenitors. Model R250 shows 
a weak and then strong peak in its light curve (Figure[l|, the 
initial peak powered by thermal energy and the second by 
the radioactive decay of ^^Ni. The heating from radioactive 
decay delays the inward-propagating recombination wave 
from ejecta cooling, regulating the electron scattering opac- 
ity (and thus the release of thermal energy), and causing the 
second peak to rise at 200-300 days, which reaches a spec- 
tacular brightness of ~ —21.5 mag. However, model R150 
produces very little ^''Ni, and therefore lacks a prominent 
second peak; the light curve is essentially thermally pow- 
ered and reaches a brightness less than that of a Type la 
SN. 

The helium core models are more compact and hence 
lack an initial thermal peak (Figure [Sjl. Model Hel30 
reaches an exceptional peak brightness of 2 x 10'*'' ergs s~', 
whereas Model HeSO demonstrates that despite being mas- 
sive and energetic, not all PISNe are bright. This steep mass- 
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200 250 300 

days since explosion 

Figure 1. Rest frame R-band light curves for the red supergiant 
progenitor models. In some models, a brief spike in luminosity 
occurs at the end plateau when radiation is released by hydrogen 
recombination. The sharpness of the spike may be exaggerated by 
the lack of numerical convergence of the ionization front recession. 




100 150 

days since explosion 

Figure 2. Rest frame R-band light curves for the helium core 
progenitor models. 



luminosity relation for PISNe suggests that to increase the 
sheer number of SNe detected, it is b etter to conduct a 
wide rather than deep survey of the sky (jWeinmann fc Lilly! 
l2005h . 

The spectra of a PISN resemble that of average SNe, 
with P-Cygni line profiles on top of a blackbody, see Figures 
O m For the RSG models, at early times, the spectrum is 
rather featureless with only weak Balmer and calcium lines, 
reflecting the low abundance of metals in unburned ejecta. 
The spectral energy distributions of the models are blue at 
earlier times (< 50 days) but become redder over time as 
the expanding ejecta cools. In addition, line blanketing of 
the bluer wavelengths becomes more prominent over time, 
as the photosphere recedes into the deepest layers which are 
abundant in freshly synthesized iron group elements. For 
PISNe at the redshifts of reionization, JWST will mostly be 
observing in the rest frame UV, so it is important to use 




2000 4000 6000 
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Figure 3. Time evolution of the rest frame spectra for the R250 
red supergiant model. 
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Figure 4. Time evolution of the rest frame spectra for the Hel30 
helium core model. 



more acc urate spectral models, rat her than the blackbody 
models of ^ Scannapieco et al.l (|2005l ). 

Spectroscopic or rest frame UV observations of PISNe 
may be able to constrain the metallicity of the progenitor 
star. However, the hydrogen envelope may be polluted by 
newly synthesized metals mixed out during the explosion. 
[Ch en et al.l (|201ll ) simulated multi-dimensional models of 
PISNe to predict the degree of mixing. They found rela- 
tively small fluid instabilities generated from burning at the 
boundaries of the oxygen shell, and concluded that PISNe 
keep their onion-shell structure in the explosion, until the 
reverse shock passes which generates Rayleigh- Taylor insta- 
bilities. This is in contrast with CCSNe, in which a shock 
runs through the inner metal-rich core, inducing the growth 
of instabilities and mixing. Also, ordinary Pop II/I CCSNe 
have non-zero metallicity in their hydrogen envelopes to 
begin with. Hence, metal lines in early-time spectroscopy 
might be able to distinguish PISNe from CCSNe, before the 
photosphere has receded deep into the ejecta. With a little 
mixing, N and possibly some C and O might appear in the 
early spectra of PISNe, but PISNe will not have any Si, Ni, 
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Fe lines (|joggerst fc WhalenllioTll ). This best applies to the 
red supergiant models, as the Helium core models undergo 
significant burning and have spectra that show many metal 
lines at maximum light. 



3 SUPERNOVAE DURING REIONIZATION 

Observations of quasar absorption spectra l|Fan et al.ll2006l ) 
indicate that reionization was completed by z = 6. It is be- 
lieved that most of the ionizing photons came from stars 
l|Loebll2010l : iBouwens et al.|[201ld ). Although the very first 
stars could have ignited as early as z ~ 30 — 40, due to the 
exponential nature of structure formation, most of ionizing 
photons originated from stars born in the later stages of 
reionization at z ~ 10. Although it is not known which pop- 
ulation of stars dominated at this epoch, an unusual stellar 
mass function 10-20 times more efficient than the standard 
Salpet er IMF in producing ionizing photons is required at 
z ~ 6 l|Cenll201ol ). This favors the existence of a top-heavy 
Pop III stars at these redshifts, which may be observable 
via their extraordinary deaths as PISNe if the IMF included 
mostly stars between 140 and 260 Mq. Moreover, the ob- 
served rates of PISNe, CCSNe, and Type la SNe may be 
used to infer the IMF responsible for reionization at z > 6. 

3.1 The Initial Mass Function 

The ionizing photon yield per baryon incorporated into 
present day stars with a Salpeter IMF is rf^ ~ 4000. How- 
ever, if the IMF is dominated by massive metal free stars 

1 M > lOOMp) ), then 97 can be up to a factor of 20 higher 
Bromm et al.ll200ll : lR"aiter et al.ll20ld '). The transition from 
Pop HI to Pop II/I star formation is thought to occur at a 
critical metallicity of Zcrit ~ 5 x 1O~*Z0, above which cool- 
ing and fragmentation become effici ent, which stops the pre f- 
erential formation of massive stars l| Bromm fc Loebll2003l ). 

This transition can be associated with the assembly of 
atomic H cool ing halos with virial temperatures > lO^K 
l|HaimarJl2009l '). Molecular hydrogen is fragile to photodis- 
sociation, and the molecular coolant in halos are likely de- 
pleted after a single episode of metal free star formation. 
Therefore, molecular hydrogen halos are unlikely to allow 
continued formation of stars above Zcrit- Subsequent star 
formation only occurs when the deeper gravitational poten- 
tial wells of atomic H cooling halos are assembled, corre- 
sponding to a virial temperature of Tvir ~ lO'^K and a 
minimum halo mass of M^aio ~ 10* Mq The gas in these 
halos will thus have already gone through a burst of primor- 
dial star formation, and contain traces of metals, leading to 
Pop II star formation. Most of the photons responsible for 
reionization will come from the profusion of these Pop II 
stars in this scenario, although without contribution from 
Pop HI stars, this may require an unrealistic star formation 
efficiency, see Figure [S] 

There is another possibility. Most molecular H2 cooling 
halos may not have formed stars at all, due to global H2 pho- 
todissociation by an early cosmic background of 11.2-13.6 eV 
photons (the Lyman- Werner band), to which the universe 
is otherwise transparent. In this scenario, the majority of 
primordial star formation will appear in atomic H cooling 
halos with Mhaio ~ 10* M©. During blowouts from repeated 



Table 2. Model parameters of the different IMFs. Here a is the 
slope of the stellar mass function, i.e. dN/dlogM oc M" , and the 
slope of the Pop II IMF flattens at M < 0.5Mq. All mass units 
are in Mq . 



IMF model 4>(M) 




upper 


a 




Pop III Salpeter 


1 


500 


-1.35 


28683 


Pop III Flat 


1 


500 





77087 


Pop II 


0.1 


125 


-1.35 


5761 



SN explosions, these halos allow most of their self-generated 
metals to be accelerated into the IGM as SN ejecta, but, in 
contrast to smaller molecular H2 cooli ng halos, these halos 
hold on to most of their interstellar gas l|Mac Low fc Ferraral 
[l99i), and can have significant Pop HI star formation. Cou- 
pled with the high ionizing efficiency of massive metal free 
stars, in this scenario Pop HI stars will make a significant 
contribution to reionization. 

Hence, we consider two mutually exclusive scenarios for 
reionization, where either Pop HI or Pop II stars reionize 
the universe; the actual star formation history of reioniza- 
tion will be in between these limiting cases. As for the IMF 
in each scenario, for massive, metal free Pop HI stars, we use 
either a Salpeter IMF slope dN/dlogM oc M"^'^'', or a flat 
IMF slope dN/dlogM oc M° hinte d by recent simulations 
l|Clark et al.ll201ll : lGreif et al.ll201ll ). with M^pp^r = 500Mq 
and Miower = IMq. Note that the resulting PISN rates are 
not sensitive to the upper and lower mass bounds of reason- 
able Pop HI IMFs. As long as Mnpper > 26OM0 , due to the 
steepness of dN/dM, there are not enough stars at the most 
massive end to affect the overall normalization. Moreover, 
for our SFR model described in the next section, any rea- 
sonable Miower ranging from 0.1 - IOMq makes negligible 
difference to the PISN rates of Pop HI stars. This is because 
we normalize star formation by requiring the generation of 
enough stellar UV photons necessary to reionize the uni- 
verse, and massive stars M > IOMq are drastically more 
efficient at producing ionizing photons. In essence, we fix 
the number of massive stars produced in any IMF, but are 
free to vary the number of low mass stars, as the latter do 
not contribute to reionization anyway. Miower ~ IMq was 
chosen to mat ch the smallest Po p HI stars seen in recent 
simulations bv lClark et al] |2oi3). 

For Population II stars forming with traces of met- 
als, we use a Salpeter IMF with Mupper = 125Mq and 
Miower = O.IMq, where we include a factor 0.7 in the mass 
integral to account for the reduced numb er of low mass stars 
in a realistic IMF (|Fukugita et al.lll998l ). compared to the 
original Salpeter IMF. The different IMF models are tabu- 
lated in Table [H 



3.2 The Star Formation Rate 

We calibrate the SFR by requiring enough UV photons are 
produced by stars so as to ionize the intergalactic medium 
(IGM) by the end of reionization. This requires C ~ 10 
ionizing photo ns per baryon in t he IGM, accounting for re- 
combinations (|Trac fc Cenll2007l ). Using the time-averaged 
ionizing flux and stellar lifetime for individual stars from 
ISchaereil(|2002il . we find the number of ionizing photons per 
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baryon incorporated into stars rj^ (M) as a function of stellar 
mass, for Pop III stars and early Pop II stars {Z — 1/50Zq). 
For a given stellar track, the average ionizing photon per 
baryon in star is thus: 



jT]-y{M)(l>{M)MdM 
/ (j>{M)MdM ' 



(2) 



where (j>{M) denotes the IMF. Then, the fraction of total 
baryons in the universe that are in stars Fs{z) must satisfy: 



Fs{Zend)'n^fe: 

c 



= 1, 



(3) 



where z^nd is the redshift at the end of Reionization, chosen 
to be Zend = 6 iu our model, and where f^sc is the escape 
fraction of ionizing photons from the host galaxy into the 
IGM. In the calibration of the SFR as a function of red- 
shift, we assume the stars instantaneously produce all the 
ionizing photons at birth that they would normally produce 
during their lifetimes. For a fixed redshift of reionization, 
this will underestimate the SFR. Nevertheless, since most of 
the ionizing radiation was dominated by the massive stars 
(M ^ IOM0), with lifetimes < 10 Myr, this is an adequate 
approximat ion . 

The mass in stars per comoving volume as a func- 
tion of redshift, p*{z), can be related to the fraction of 
gas in halos which converts to stars, i.e. the star forma- 
tion efhciency using th e Sheth-Tormen mass function 

(|Sheth fc Tormenlll999l ) of halos: 



O. 
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Pop II 
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Pop III Flat 
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Figure 5. The SFR for Pop II or Pop III stars required to 
reionize the universe by z^^d = 6. The calibrated star forma- 
tion rate is /* = 0.3%, 0.8%, and the unrcalistically high 24.1% 
for the Pop III Flat, Pop III Saipeter, and Pop II IMF models, 
respectively, suggesting that Pop II stars could not have driven 
reionization by themselves. The SFR inferred by I Bouwens et al] 
(2011b) from integrating the observed galaxy UV luminosity den- 
sities to Mj\B 18 is plotted for comparison; their substantially 
lower SFR is not surprising, as the contribution from the very 
steep faint-end slope of lower-luminosity galaxies was omitted 
llBouwens et al.l2011al V The SFR for our models and the resulting 
SN rates all linearly scale with C and /^J,. C = 10 and fesc = 0.1 
were used throughout this paper. 



p^{z) = Fs(z)pb = /, 



dM 



(4) 



Here pb is the cosmological baryon density, and Mmin ~ 
10^ Mq for both Pop II and Pop III scenarios, corresponding 
to halos with atomic H cooling. For cosmological parame- 
ters used in generating the Sheth-Tormen mass function, we 
adopt h = 0.71, r2„ = 0.27, f^A = 0.73, and fit = 0.045, 
where h is the Hubble constant in units of 100 km s~^ and 
Qm, f^A, and ilf, are the total matter, vac uum, and bary- 
onic densities in units of the critical density (jKomatsu et al.l 
I2OIII ). Since Fsiz^nd) is known via equation by plug- 
ging z = Zend iuto cquatiou Q, we can calibrate the value 
of /*, and evaluate p,{z) at any redshift. Although /, will 
generally vary with redshift, here we take /, as a constant 
for simplicity of calibration. The star formation rate is then 
simply: 



SFR{z) = 



dpf[z 
dz 



(5) 



Figure [5] shows the resulting star formation rates, using 
C = 10 and fesc = 0.1. Observations of Lyman- break galax- 
ies around z ~ 3 suggests that fesc could be larger at 
higher redshifts (jSteidel et al.|[200l] ). However, theoretically 
the high gas densities at the redshifts of the first galax- 
ies could decrease th e escape fraction down to fesc 0-01 
l|Wood fc Loebll2000l ). in which case using fesc = 0.1 is a 
conservative choice that may underestimate the SFR and 
the corresponding SN rates. To consider different choices of 
these parameters, note that the SFR in our model linearly 
scales with C and fe. 
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Figure 6. The rate R{z) of PISNe, CCSNe, and Type la SNe 
for our different IMF models for the stellar population responsible 
for reionization, per year per Mpc^. 



3.3 Snapshot Rate with JWST 

The James Webb Space Telescop^ (JWST) will include a 
Near Infrared Camera (NIRCam), with a spectral coverage 
from 0.6-5 /im with ~10 njy sensitivities in lO** s of inte- 
gration time (IOct); a Near Infrared Spectrograph (NIRSpec) 
which operates at approximately the same wavelength range. 

http://www.jwst.nasa.gov/ 
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The Mid InfraRed Instrument (MIRI) covers 5-27 /im, but 
is an order of magnitude less sensitive than NIRCam. Since 
isolated Pop III stars are likely beyond the reach of JWST, 
to test the prediction that metal-free stars had a top-heavy 
IMF ( which has been recently debated, see lHosokawa et al.l 
lj201l| )). we can either observe the cumula tive properties of 
the fi rst stars by imaging Pop III galaxies feackris son et all 
or detect their deaths as extraordinary bright super- 
nova. 

The number of new events at a given redshift that can 
be observed per unit solid angle is jWoods fc Loebl flQQS'): 

N{z) = R{z) (1 + z)~'^ for Z < Z^ax{F;v), (6) 

where Zmax{F;v) is the maximum redshift at which a source 
will appear brighter than limiting flux F at an observed fre- 
quency V, R{z) is the event rate per unit comoving volume, 
and dVc is the cosmology-dependent comoving volume ele- 
ment corresponding to a redshift interval dz. The above ex- 
pression includes the (1-l-z) reduction in apparent rate owing 
to cosmic time dilation. 

The 'snapshot rate', i.e. the total number of events (not 
per unit time) observed at limiting flux F is: 

f'°° dV 
N{F-u) = J^ dzR{z)t{z-F-u)-^, (7) 

where t{z; F; u) is the rest-frame duration over which an 
event will be brighter than the limiting flux F at redshift z 
for an observed frequency u. We find this duration from our 
spectral time series calculated with SEDONA. There is an 
implicit (l + z) factor in equation ([T]) due to the time dilation 
of the light curve, but that cancels with the {l + z)~^ reduc- 
tion in apparent rate. Although t[z; F; v) of PISNe will gen- 
erally be longer for more massive progenitors, the snapshot 
rate is not necessarily dominated by the highest mass stars, 
as they are less numerous, see Figure[7l For CCSNe, it is not 
clear how the brightness of a Type IIP SN should depend 
on the mass of the progenitor star; here we use light curves 
of a ISMp) re d giant progenitor generated with SEDONA 
asen fc Wooslev..2009il . whose broadband light curves and 
spectra agree very well with observed Type IIP SNe, which 
are observed to be the most common, at least in the nearby 
universe. However, this single CCSN model means we do 
not capture the variation in CCSN peak flux from different 
progenitors, which we do so for PISNe. 

PISNe and CCSNe occur for stars with main-sequ ence 
masses betweenj~14O-26OM0 (|Heger fc Wooslevlbooj ) and 
8-25M0 (|Smartt|[2009l ). respectively. The fate of stars be- 
tween 25-14OM0 is uncertain; due to fallback, progenitors 
more massive than ~ 4 0M(;) may fo rm black holes directly 
with no SN explosion l|Frveij [l999l ). Notably, stars in the 
mass range 95-13OAf0 may reach the pair production in- 
stability in the core, but the thermonuclear explo sion is in- 
sufficient to unbind the star (|Wooslev et al.ll2007l ). and the 
star undergoes pulsations of matter ejecta which may pro- 
duce a very bright light curve when the shells of ejected 
matter collide with each other, before the star dies as a nor- 
mal CCSN. The resulting pulsation pair-instability super- 
nova can be ultra-luminous and are presumably detectable 
by JWST. However, we do not consider such events here. 

Using a progenitor mass range of 8-25M0 to calculate 
the CCSN rates is likely an underestimate; for detailed pre- 
dictions on the number of CCSNe detectable by JWST at 



Table 3. Snapshot rate in NIRCam's field of view (10 arcsec^) on 
JWST, using the same survey parameters as Figure[7] This is the 
total number of PISNe and CCSNe in each NIRCam snapshot, 
integrated across 2^6, for the different IMF models. One can 
multiply the values below by 360 to get the snapshot rate per 
deg^. It is not clear what fraction of PISNe explode from red 
supergiants (RSG) versus exposed helium cores; the rates shown 
in the RSG and He Core columns assume all PISN explode as 
that type. Here the snapshot rate of PISNe from red supergiants 
is higher than the helium core model due to the longer duration 
of the former. The snapshot rate for CCSNe in our Pop II IMF 
scenario is in good a greement with the hig h end estimate of 24 
SNe per field found by|M esinger et alj l|2006h under similar survey 
parameters. 



IMF model 


He Core 


RSG 


CC 


Pop HI Salpeter 


0.28 


0.42 


10.43 


Pop HI Flat 


0.74 


1.03 


0.64 


Pop II 








31.83 



Table 4. Snapshot rate in MIRFs field of view (2.35 arcsec^) on 
JWST, using the F770W filter (77000 A) with 3 x lO^s exposure 
(5(t). The resulting snapshot rate is an order of magnitude worse 
than NIRCam; however, these results suggest that MIRI can be 
used as a follow-up instrument to distinguish bright PISN events 
from core-collapse events. Since MIRI is much less sensitive, the 
brighter He core models are more readily observable, while CC- 
SNe cannot be seen at all in this integration time. 



IMF model 


He Core 


RSG 


CC 


Pop HI Salpeter 


0.05 


0.04 





Pop HI Flat 


0.15 


0.12 





Pop II 












the redshifts of reionization, see lMesinger et all (|2006h . who 
also take into account the variation in peak magnitude of 
observed CCSNe and the effects of dust extinction. For a 
fixed progenitor mass range, we calculate the SN rate per 
comoving volume R[z) and find the snapshot rate shown in 
Tables [S] [4] and Figures [6] [T] In Figure [71 we use each red 
supergiant model as a proxy for the light curves of all pro- 
genitors similar in mass, e.g. R250 represents all progenitors 
in mass range 226-26OM0. 

If both Pop HI and Pop II stars contributed to reioniza- 
tion, the actual IMF will be a mixture of the Pop HI and Pop 
II IMF used above. By counting the number of each type of 
SN found in JWST snapshots, the IMF of these early stellar 
populations can be constrained, and the relative contribu- 
tion of Pop HI and Pop II stars toward reionizing the uni- 
verse can be inferred. To reduce selection effects due to the 
different intrinsic luminosity of the SN, a sufficiently deep 
exposure is needed, to enable observations of both types of 
SNe at the peak of their light curves should they exist at 
the target redshift. 

For the integration time assumed in Figure [T] one can 
directly characterize the ratio of PISNe to CCSNe before 
2 ~ 8, and set existence limits on top heavy Pop HI stars 
up to z ~ 10 with a survey of ~ 10 JWST fields. 
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6 7 8 9 10 11 12 13 14 15 6 7 8 9 10 11 12 13 14 15 

Redshift Redshift 

(a) Pop III Salpeter IMF model (b) Pop III Flat IMF model 

Figure 7. Differential snapshot rate in the 10 arcsec^ field-of-view of NIRCam on JWST, caleulated using the sensitivities of 

the F444W filter (44400 A) with t = 3 X 10*s integration time at 5cr (a flux threshold of 2 njy). At this sensitivity, each type of SN 
appears in rough proportion to their actual event rates (Figure [B]! up to z ~ 8, past which the ratio of detected PISNe versus CCSNe 
turns over, with only the brighter PISNe staying in view. As we have not accounted for the intrinsic scatter in the luminosity of CCSNe, 
the actual turnover will be less sharp. A similar turnover exists for the less massive PISN models R175, R150 starting at z 11, past 
which the more massive progenitors are more likely to be seen in the field of view despite be ing less numerous in number given the IMF. 
The snapshot rate of CCSNe in Figure [Tal is less than that implied in lMesinger et all | |200(J) , as their assumed IMF is closer to our Pop 
II IMF model. 



3.4 Probing Intermediate Mass Stars with Type 
la SNe 

To probe the intermediate mass range (~ 1 — 8Mq) of the 
IMF during reionization, one may use Type la SN rates. 
Type la supernovae (SNe la) are thought to occur when a 
white dwarf nears the Chandrasekhar mass, resuhing in a 
thermonuclear explosion. This requires the white dwarf to 
accrete mass from a binary companion. Although the physics 
behind SNe la have been widely studied using both observa- 
tions and theoretical simulations, there is still no consensus 
on the mechanisms that proceed the supernova. The single 
degenerate model proposes the companion to be a main se- 
quence or giant star, which donates mass via Roche lobe 
overflow, whereas the double degenerate scenario considers 
the merger of two white dwarf stars; the latter may be nec- 
essary for at least s ome observed Type la SNe dSloom et al.l 
I2OI2I : ISchaefer fc Pagnotta ,2012. ). Either way, after stellar 
birth it takes the main sequence lifetime of the progenitor 
star plus an additional delay time for the Type la SN to 
proceed. 

Hence, the rate of SNe la is empirically parametrized 
to follow the star formation rate (SFR), but shifted toward 
lower redshift after taking the delay time into account. The 
SN rate at a redshift z or cosmic time t, R{z) — R{t), is 
given by a convolution of the SFR over delay times, 

R{t) = / SFR(t - r)DTD(r)dr, (8) 
Jo 

where DTD(r) is the delay time distribution (SNe per unit 
time per unit stellar mass formed), in which r is the time 
elapsed between the formation of the progenitor star and 
the explosion of the SN la. Note that since the DTD(t) is 
normalized to the total stellar mass formed, it only indirectly 
reflects the physical efflciency of SNe la fr om their actual 
progenitors of 3 — 8M@ stars (jNomoto et al.lll994i '). 



In previous reionization literature (iHaimanl bOO^l.Tvpe 
la SNe were expected to be extremely rare at high redshifts 
{z > 6), as the delay between the formation of the progenitor 
and the SN event was thought to be longer than the age of 
the Universe at these redshifts. However, this view should 
be reconsidered in light of recent converging eviden ce for 
a prompt population of S Ne la, see recent work by iMaoj 
(j2010l ). [Graur et all ()201ll ). and references within. 

IScannapieco fc BildstenI (|2005l ) calibrated the prompt 
rate via the 'B' parameter, a constant of proportionality be- 
tween the SFR and the prompt SN la rate, equivalent to 
the number of prompt SNe per unit stellar mass formed. 
The delayed component is characterized via the parameter 
'A' which is the constant of proportionality between galaxy 
mass and the delayed SN la rate. We ignore the A compo- 
nent as this delay exceeds the age of the universe during 
reionization. The value of B is calibrated at low redshifts, 
for example B = 2.7 - 11 x 10~^M~^ in iMaoz fc Badenej 
(2010), for prompt delay times T £ (35,330) Myr; here 
we adopt B = 3 x 10~"^A/g^. From B, we set a uniform 
DTD(r) = B/AT. Using the Pop II SFR and this DTD(r), 
we calculate the event rate R{t) of Type la SNe shown in 
Figures [6] and [H Since the validity of these estimates de- 
pend on DTD(r), we assume that the astrophysics involved 
in shaping the forming efficiency and delay time of Type la 
SNe is not very sensitive to the cosmological epoch. 

As seen in Figure [8l in the scenario where Pop II stars 
dominated reionization, we expect ~1 new SNe la every year 
per unit redshift at z ~ 6 — 7 in NIRCam's field of view of 
~ 10 square arcseconds. As the AB magn itude of SNe la 
during peak ( Hillcbrand t fc Nieme vcr 2000) is Mb ~ Mv ~ 
— 19.30 with a dispersion of 0.3, several hours integration 
time on JWST will be sufficient to catch a Type la SN near 
peak flux at these redshifts. 

Regardless of which population of stars dominated 
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Figure 8. Type la SNe rate in NIRCam's field of view, per year 
of observation per dz. For the Pop II IMF model for reionization, 
~1 new SNe la will occur every year per unit redshift at z ~ 6 — 7 
in NIRCam's field of view. The 20 days around the SN la peak 
flux at that redshift is equivalent to half a year in observer frame. 
Every NIRCam snapshot of the sky, with sufficient integration 
time (e.g. 3 X 10** s to get absolute magnitude above -18 at z = 8 
with signal-to- noise ratio of 5), will have a ~ 50% probability of 
finding a Type la SN near peak flux. 



reionization, Type la SNe offers a way to probe the interme- 
diate mass range of the reionization IMF. In addition, the 
z > 6 sky offers an unambiguous way of isolating a prompt 
population of SNe la, as the universe was not sufficiently 
old in the epoch of reionization for the delayed component 
of Type la SNe to contribute any events. Therefore, Type la 
SNe at 2 > 6 could be used to test whether the delay times 
are indeed connected to SN la formation mechanisms and 
properties. Finally, the existence of standard candles during 
the epoch of reionization could be useful for cosmological 
measurements. In the Appendix, we discuss that probing 
reionization history with Type la SNe by constraining the 
global ionization fraction using Thomson optical depth mea- 
surements requires an unrealistic survey time for JWST. 



3.5 Typing the Supernovae 

At lower redshifts, SNe are usually typed by spectral lines. 
Furthermore, the smoking gun evidence for a PISN is the 
measurement of a large core (> 50 Mq) composed of he- 
lium or other heavier elements. In the case of SN 2007bi, 
iGal-Yam et"aLl l|2009l ) analyzed the nebular spectrum 16 
months after peak light to infer ~ 4Mq of ^^Ni, implying 
a large core mass > SOMq as i n a pair-instability explo- 



sion (but see lMoriva et al.l (I2OI0I ) for a CCSN model for SN 
2007bi that ejects 6.IM0 of ^^Ni). However, spectroscopic 
typing of high redshift SNe seen by JWST may be unreal- 
istic; for example, at around 3 jj,m on NIRSpec, achieving a 
signal-to-noise ratio of 5 at redshift ~ 8 would require one 
full day of integration time. Alternatively, in anticipation of 
Pan-STARR^and LSStEI increasing the number of photo- 



^ http:/ /pan-starrs. ifa.hawaii.edu/public 
^ http://www.lsst.org/lsst/ 



metrically detected SNe to a few hundred thousand in the 
next two decades, much work has been done in the photo- 
metr ic identification and classification of SNe (jKessler et al.l 
I2OI0I ). As photometric classification of SNe matures, it could 
complement or replace the spectroscopic typing of high-z 
SNe, reducing the required JWST time. 

The hght curves of the more massive PISNe are very 
luminous (10*^-10'*'* erg s~^) and long-lasting (~ 300 days), 
characteristics that do not exist for most other types of SNe. 
As long as the SN redshift is known, multi-epoch observa- 
tions can determine its rest-frame luminosity and duration, 
and identify the more massive PISNe explosions. Aside from 
the most energetic events, typing PISNe using their magni- 
tude and color will be difficult. Despite their enormous ki- 
netic energies of ~ 10^^ ergs, the peak optical luminosities of 
PISNe are similar to those of other SNe, even falling below 
the la and II curves for smaller mass progenitors. The ma- 
jority of PISNe will actually be these dimmer events. Also, 
since PISNe spend most of their lives in the same tempera- 
ture range as other SNe, their colors are also similar. 

An extended light curve, rather than an extreme lumi- 
nosity or unusual color, may therefore be the most important 
signature of PISNe. In particular, the distinguishing feature 
of PISNe is its exceptionally long rise time, S> 100 days in 
the rest frame. Also, the detection of a slow decline rate that 
follows the decay rate (~0.01 mag/day) of ^®Co, the prod- 
uct of ^®Ni decay, would provide strong evidence that the 
SN synthesized significant amounts of ^^Ni. At z = 8, even 
with time dilation, this results in ~0.4 mag variation per 
year, which should be within the sensitivity of a multi-year 
JWST survey; the decline in the bluer bands are 2-3 times 
larger, due to the onset of iron group line blanketing. The 
detection of a secondary maximum in the light curve also 
supports the synthesis of ^''Ni. However, the lack of a sec- 
ondary peak does not rule out a large presence of ^ ^Ni, as 
stron g radial mixing could smear out the two bumps (jKasenl 
l2006h . 



3.6 Survey Strategies 

The long duration of high redshift PISN light curves, pro- 
longed by cosmological time dilation, poses a great challenge 
for detecting them as transients. At z~7 the light curve of 
a PISN can last for over 1000 days in the observer's frame. 
Without spectroscopic measurements, the telltale sign of a 
massive progenitor PISN is an incredibly long plateau in its 
light curve. Therefore, instead of a threshold experiment, 
we suggest a search strategy that involves taking a series 
of 'snapshots' of a field, each snapshot separated by ~ 1 
year, and searching for variations in the flux of objects in 
successive images. Since Pop HI star formation occurs in the 
smallest galaxies, blank-field surveys should be the sufficient 
for searching for PISNe. 

The total number of SNe detected in a survey of total 
integration time tsurv is 

1 tsurv ^^FOV 

2 tea 



47r 



(9) 



where Aiipov is the instrument's field of view, tsurv/texp 
is the number of fields which can be tiled within the survey 
time tsurv, and N^xp is the snapshot rate from equation (O, 
i.e. the number of SNe bright enough to be detected in an 
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exposure of duration texp l|Haimanll2009l ). The factor of i 
is included to account for observations in 4 color bands (2 
pairs of filters, as NIRCam observes in two bands simultane- 
ously using a dichroic) for determining photometric redshift 
and typing of the SNe. To detect SNe by their variability, 
each field requires repeated observations, and therefore any 
survey should piggyback on fields that have already been 
observed. 

In the case where several fields are already available 
from other JWST surveys, the snapshot rates given by Ta- 
ble [3] suggest that a dedicated, long program may not be 
required to detect dozens of high redshift SNe. Due to the 
order-of-magnitude difference in the snapshot rate of PISNe 
vs CCSNe for the different IMF models, more than 10 fields 
with foUowup repeated imaging should already help con- 
strain the stellar population responsible for reionization. 
Cosmic variance will affect the total number of SNe for small 
number of fields, but the ratio of PISNe to CCSNe would still 
be indicative of the IMF. For example, if Pop III (Flat IMF 
model) and Pop II stars had equal contribution to reion- 
ization (which means Pop II stars dominate Pop III stars 
by roughly 20-to-l in total mass), one could use 20 images 
conducted for other programs as references and only revisit 
the same image twice for a total of 3 snapshots per field of 
view. Observing in 4 bands, for a total of 28 days integra- 
tion time over 2 years, such a survey expects to see ~10 red 
supergiant PISNe and ~300 CCSNe. 

To the extent that PISN spectra can be represented as 
a distribution of blackbodies at different temperatures, since 
the temperature and redshift would be degenerate, it will be 
impossible to acquire photometric redshifts without further 
information about the SN epoch. However, our simulated 
spectra show significant deviations from a blackbody in the 
UV {I < 3500A) due to metal-line blanketing in the SN 
photosphere, providing spectral and photometric signatures 
that could be used as redshift indicators, depending on their 
strength. 

3.7 Luminosity Function 

Although the UV flux of PISNe is relatively short lived, the 
more massive PISNe stay bright in its rest frame visible band 
for over a year. Given this brightness and long intrinsic dura- 
tion, coupled with the (l + z) time dilation at high redshifts, 
it is conceivable that PISNe could contribute to the lumi- 
nosity function of all objects at high redshifts when galax- 
ies were dim. Figure [9] illustrates the luminosity function of 
PISNe at ~4000A, calculated using the helium core progeni- 
tor models for PISN luminosity, and the Pop III Flat or Pop 
III Salpeter models for the star formation rate. Shown for 
comparison are the projected galaxy luminosity functions 
at high redshifts, using the Bouwcns et al. (2Cnibj) best fit 
Schechter parameterization for the UV luminosity function, 
and shifting to the visible band using U — V ~ 0.4, 0.3 for 
z = 7, S respec tively, measured using the Spitzer Infrared 
Array Camera l|Labbe et all l2010al lbl). Applying this U-V 
shift is a crude approximation, as luminous and faint galax- 
ies have different rest frame UV-to-optical color; however, 
we are most interested in the bright end of the luminosity 
function, where this current U-V measurement is applicable. 

The luminosity function for PISNe implied by our Pop 
III IMF models overlaps with the galaxy luminosity function 




Figure 9. Luminosity functions of PISNe at 4000Afor 2 = 7 
and 8, for the Pop III Flat IMF or Salpeter IMF star formation 
rate models, compared with the galaxy luminosity function. For 
the Pop III Flat model, the volumetric count of PISNe exceeds 
galaxies past M^g ~ —21.5; the dominance of PISNe should 
become greater at higher redshifts, as galaxies decrease in lumi- 
nosity while PISNe stay the same. The PISN luminosity functions 
at 2 = 7, 8 overlaps coincidentally because the increase in time 
dilation compensates for the decrease in PISN event rate. 

at the brightest magnitudes. If a top-heavy Pop III IMF was 
solely responsible for reionization, PISNe will contaminate 
the brightest end of the galaxy luminosity function, unless 
great care is taken to remove these supernovae. Since the 
volumetric count of the brightest galaxies and PISNe is very 
low, it will take a wide infrared survey to observe this effect. 



4 DISCUSSION 

In our discussion we ignored complicating factors such as 
metallicity and rotation, and calculated the PISN and CCSN 
event rate using only the SN progenitor mass range along 
with the star formation rate. However, at low redshifts z < 1, 
the measured CCSN rate is a factor of ~ 2 smaller than that 
predicted by the analogous calculation using the measured 
cosmic star formation rate. The discrepancy is likely due 
to many intrinsica lly low-luminosity or obscured SNe being 
missed in surveys (|Horiuchi et al.|[201ll ). As this discrepancy 
is lower than the uncertainty in our SFR model parameters, 
and we already account for lower intrinsic luminosities for 
the lower progenitor mass PISNe, we do not take obscuration 
into account for our predictions of the SN rate as seen by 
JWST. 

The IMF of early stellar populations responsible for 
reionization should also leave an imprint on the metal en- 
richment pattern via their SN products. So far, the abun- 
dance patterns observed t o date in extremely meta l-deficient 
stars in the Galactic halo (|Beers fc Christhebll2005l ) are more 
consistent with an IMF that produced much more CCSNe 
instead of PISNe (iJoggerst et al.ll201Gl ). However, in previ- 
ous surveys, subtle selection effects might have disfavored 
finding PISN-enriched stars; the metal yields of PISNe are 
so high that the metal a bundances of stars formed out of 
PISN ejecta l|Greif et al.li2008i ') are already higher than the 
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metallicity ranee tar geted by metal-deficient star surveys 
l|Karlsson et al.li2008l ). 



Large carbon eniiancements observed in metal-poor 
stars, when interpreted as the outcome of pollution by 
winds from binary companions that have gone through 
the AGB phase, suggest the existence of a large num- 
ber of intermediate-mass stars (~ 1 — 8Mq) at high red- 
shifts (Tumlinson ^007a,b). Alternatively, nucleosynthesis 
in faint CCSNe from higher mass stars could also ex- 
plain the observed ca rbon enhancement in metal-poor stars 
l|lwamoto et al.|[2005l ). Observing the Type la SN rate dur- 
ing the epoch of reionization will be an complementary 
way to test these models, and constrain the number of 
intermediate-mass stars at high redshifts. 

The predicted initial mass range of ~140 to 260 Mq 
for PISN progenitors as sumed the stars to be non-rotating 
jHeger fc Wooslevll2002l ). However, observat ions find that at 
very low metallicities, stars rotate faster (|Martavan et al.l 
I2OO7I ). The fast rotation of the first stars is supported 
by the la test hydrodynamic simulations of their formation 
I Stacy et al. 2011), and also by observations of anomalously 
high abundanc es of Ba and La with r espect to Fe in ancient 
low-mass stars jChiappini et al .l201ll '). which could originate 
in metal-poor fast-rotating massive stars. Generally, rota- 
tion should increase the required P ISN progenitor mass by 
increasing mass loss. iMevnet et al. found that, con- 

trary to the usual M oc Z^'^ scaling relation, rotating stars 
at very low metallicity Z ~ 10~^ to 1O~*Z0 show a large 
mass loss, up to ~ 50%, mainly resulting from efficient mix- 
ing of stellar nucleosynthesis produ cts into the stellar sur- 
face. However, lEkstrom et al.l (|2008t ) found that for strictly 
Z = stars, the mass loss is very low, even for models that 
reach critical velocity in the main sequence. These results 
imply that, for rapidly rotating Pop HI stars to die as a 
PISN, the required progenitor mass is extremely sensitive 
to whether the star is truly metal- free or not. 

At much lower redshifts, PISNe have likely already been 
observed, most persuasively in the case of the very lumi- 
nous and long duration event SN 2007bi (|Gal-Yam et al.l 
I2OO9I ). Other more recent candidates include PTF lOnmn 
(Gal- Yam in preparation; Yaron et al. in preparation) and 
PSl-llap (Rubina Kotak et al. in preparation). As pris- 
tine gas was recently observed at redshifts after reioniza- 
tion IjFum agalli et al. 2011), it is possible that some low red- 
shift PISNe have Pop HI progenitors born out of surviving 
pockets of metal-free g as; the rates of PISNe i n this sce- 
nario was considered bv lScannapieco et al.l (|2005l ). However, 
the metallicities of the host galaxies of SN 2007bi and PTF 
lOnmn are well ab ove the metallicity threshold required to 
form Pop HI stars (| Young et al.ll2010l ). Therefore, it is plau- 
sible that PISNe can have very massive Pop II/I progenitors 
as well, perhaps born via the merger of stars in collisi on run- 
aways in young, dense star clusters (|Pan et al.ll201ll ). 



5 CONCLUSIONS 

We analyzed simulated light curves and spectra of pair- 
instability supernovae for a variety of progenitor masses and 
envelope types, and found that the supernovae from the 
more massive progenitors are super-luminous and have ex- 
tended light curves, traits that would help photometrically 



distinguish pair-instability supernovae from other types of 
supernovae using repeated snapshots. We calculated the 
rates and detectability of pair-instability, core collapse, and 
Type la supernovae during the redshifts of reionization, and 
showed that it is possible to constrain the initial mass func- 
tion of stars at that time, and identify the stellar population 
responsible for reionization. If Pop HI stars made the dom- 
inant contribution of ionizing photons during reionization, 
the bright end of the galaxy luminosity function will be con- 
taminated by pair-instability supernovae. 
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APPENDIX A: ON THE DIFFICULTY OF 
MAPPING REIONIZATION HISTORY WITH 
TYPE lA SNE 

Al Constraining the Ionization Fraction 

Barring some extreme evolution of the IMF, the neutral frac- 
tion of the IGM is expected to rise quickly toward high red- 
shift, with the mean neutral fraction of the IGM expected 
to reac h 6-12% a t z=6.5, 13-27% at z=7.7 and 22-38% at 
z=8.8 |Cenll2010l ). Currently, the most stringent observa- 
tional probe on the ionization history of the IGM is the total 
Thomson optical dept h seen by WMAP, r = 0.088 ± 0.015 
iKomatsu et al.lbOHl ). The Planck mission is projected to 
reduce the error bars to 0.01. However, r cannot break de- 
generacies between different reionization histories; for ex- 
ample, both a rapid, early reionization or an extended, 
late reionization may have the same total Thomson opti- 
cal depth. Finding r(z) using Type la supernovae at high 
redshifts would break this degeneracy. 

We set up a toy model of the global average ionization 
fraction X{z) using the Fermi-Dirac form for the ionization 
fraction: 



X{z) 



+ 1 



(Al) 



where Zre, A are model parameters that characterize the 
redshift and duration of reionization, respectively. Then 
for an luminous object at redshift Zobs, the Thomson elec- 
tron scattering optical depth is the integral of X {z)neO'T , 
the ionization fraction times the electron density times 
the Thomson cross sectio n integrated along proper length 
(|Shull fc Venkatesanlliooi ). 



r{zobs) 



X{z)ne(TT- 



-dz 



(A2) 
(A3) 



{l + z)H{z) 

+ THelll X Q{Zobs — ZHelll)- 

Here H{z) = Ho[nrr,{l + zf + flA]^^^ ■ The second term with 
THeiii — 0.002, and O as the Heaviside step function, comes 
from the full reionization of Hell to Helll around ZHeiii ~ 3 
l|Shull et al.ll2004 ). 

Using this model with hypothesized optical depth mea- 
surements, we use Bayesian methods to find the corre- 
sponding probability distribution for our reionization his- 
tory model parameters Zre , A shown in Figure lAll These 
measurements along with the known optical depth to the 
CMB can help constrain the duration of reionization A. 

A2 Survey Feasibility with JWST 

Using the Phillips relation l|Phillipsl ITooil '). we can utilize 
Type la SNe as standard candles by measuring the shape of 
the light curve. Specifically, to characterize the light curve 
and find Ami5{B), i.e., the decline in the B-magnitude light 
curve from maximum light to the magnitude 15 days after 
B-maximum, we should to take a snapshot every 3 days 
for roughly 20 days (5 days before and 15 days after peak), 
which at redshift 8 means returning to the same field of view 
once every month due to time dilation. Note that above 
z > 10, the B ban d (~ 4 400A) redshifts out of NIRCam. 
Also, iBrandt et al.l (|201G| ) found that the prompt channel 
Type la SNe are more luminous (high-stretch, slow declin- 
ing), and thus have a lower Ami5(B). For NIRCam, we can 
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Figure Al. Onc-sigma (68%) contours of Zre, A given measure- 
ments of T^Zohs = 8) ^■ud CMB total optical depth. Flat priors 
of Zre = 5 — 15, A = 0.25 — 3 are assumed for the x,y axis, re- 
spectively. The straight, long-dashed, and short-dashed line rep- 
re sent the contours fr om measurements of tcmb = 0.088 ± 0.015 
(Komatsu et alj|201ll) . T(zai,s = 8) = 0.07 ±0.01, and T(Zobs = 
8) = 0.04 ± O.Ol, respectively, the latter two of which are hy- 
pothetical values of optical depth we may measure from a large 
sample of SNe la. 



see up to -18 AB magnitude at redshift ~ 8 with 3 x 10* 
s integration time (~ 8 hours), with a signal-to-noise ratio 
of 5, and fully characterize the SN light curve as long as 
Ami5(B) < 1.3. This is roughly ~60 hours of total integra- 
tion time over half a year for each field of view with a Type la 
SN. Unrealistically assuming command of all JWST's time, 
the light curve of ~150 Type la per year can be fully mapped 
out. 

As seen in Figure HJ JWST might see ~1 Type la SNe 
at z = 6 — 7 for every snapshot it takes with NIRCam. To 
find interesting results about reionization history, one should 
probe the ionization fraction before the end of reionization, 
at z ~ 8 or above. Here, the SNe are dimmer and the rate 
is smaller, though by only a factor of ~ 2 each. A bigger 
obstacle is, at this redshift, many color bands redshift out 
of NIRCam's range (at z ~ 8, only the UBV bands at the 
SN rest frame are still accessible). Since the calibration of 
Type la SNe magnitudes relies on multiple color bands, it 
is not clear the often quoted ~ 0.20 standard deviation in 
distance modulus can be achieved. Moreover, calibrating the 
intrinsic luminosity of SNe la down to AM ~ 0.20 is equiva- 
lent to a variation in optical depth of At ~ 0.18. To get the 
Thomson optical depth to 1% precision at a fixed redshift 
bin, required to make meaningful constraints on reionization 
history, would require over 300 independent Type la SNe, 
or over two full years of JWST's integration time. This is 
clearly not feasible. 
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